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A Deficit of High-Redshift, High-Luminosity X-Ray Clusters: 
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ABSTRACT 

From the Press-Schechter mass function and the empirical X-ray cluster luminosity- 
temperature (L-T) relation, we construct an X-ray cluster luminosity function that can 
be applied to the growing number of high-redshift, X-ray cluster luminosity catalogs 
to constrain cosmological parameters. In this paper, we apply this luminosity function 
to the Einstein Medium Sensitivity Survey (EMSS) and the ROSAT Brightest Cluster 
Sample (BCS) luminosity function to constrain the value of tt m . In the case of the 
EMSS, we find a factor of 4 — 5 fewer X-ray clusters at redshifts above z = 0.4 
than below this redshift at luminosities above Lx = 7 x 10 44 erg s^ 1 (0.3 - 3.5 keV), 
which suggests that the X-ray cluster luminosity function has evolved above L*. At 
lower luminosities, this luminosity function evolves only minimally, if at all. Using 
Bayesian inference, we find that the degree of evolution at high luminosities suggests 
that O m = 0.96la32, 

given the best-fit L-T relation of Reichart, Castander, & Nichol 
(1998). When we account for the uncertainty in how the empirical L-T relation evolves 
with redshift, we find that £l m ~ 1.0 =t 0.4. However, it is unclear to what degree 
systematic effects may affect this and similarly obtained results. 

Subject headings: cosmology: observation — cosmology: theory — galaxies: clusters: 
general — galaxies: luminosity function, mass function — X-rays: galaxies 

1. Introduction 

The combination of the Press-Schechter mass function (e.g., Press & Schechter 1974; 
Oukbir & Blanchard 1992; Lacey & Cole 1993) and present and future X-ray cluster catalogs 
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presents a unique opportunity to constrain the cosmological mass density parameter, £l m . The 
Press-Schechter approach offers a number of advantages over various, more traditional methods 
of measuring this parameter. First of all, numerical simulations reproduce the Press-Schechter 
mass function to a high degree of accuracy (e.g., Eke et al. 1996; Bryan & Norman 1998; Borgani 
et al. 1998). Secondly, unlike methods that only probe Q m over small spatial scales — methods 
that may be insensitive to an underlying, more uniformly distributed component of the dark 
matter — the Press-Schechter approach probes Vt m over the scales of X-ray cluster catalogs, which 
now have limiting redshifts of about unity. Thirdly, the Press-Schechter approach appears to be 
relatively insensitive to a cosmological constant (e.g., Henry 1997; Mathiesen & Evrard 1998; Eke 
et al. 1998; Viana & Liddle 1998); consequently, a Press-Schechter-based determination of Q m 
might be compared to independent determinations of the deceleration parameter, for example, 
to constrain the cosmological constant. Finally and most importantly, a number of independent, 
high-redshift, X-ray cluster luminosity catalogs with well- understood selection functions are and 
will soon be available (see §4). We discuss potential problems with the Press-Schechter approach 
in §4. 

The archetypal X-ray cluster catalog is the X-ray cluster subsample of the Einstein Medium 
Sensitivity Survey, which we refer to here as the EMSS. A complete description of this sample 
and its selection criteria can be found in Henry et al. (1992) (see also Gioia et al. 1990b; Stocke 
et al. 1991; Gioia & Luppino 1994). The original EMSS consists of 93 X-ray clusters, of which 67 
have redshifts z > 0.14. Nichol et al. (1997) updated this z > 0.14 subsample with 21 ROSAT 
luminosities and optical information from the literature. This revised EMSS consists of 64 X-ray 
clusters, of which 25 have been updated. The redshift and luminosity ranges of the revised EMSS 
are 0.140 < z < 0.823 and 7.5 x 10 43 < L x < 2.336 x 10 45 erg s" 1 (0.3 - 3.5 keV). The EMSS is of 
great importance because at present, it is the only X-ray cluster catalog that probes masses above 
M* at high redshifts, where the Press-Schechter mass function is the most sensitive to Q. m (see §3). 

Assuming that X-ray clusters correspond to virialized, dark matter halos, the Press-Schechter 
mass function describes how the X-ray-selected cluster mass function evolves with redshift. 
Unfortunately, X-ray-selected cluster mass catalogs that span sufficiently broad ranges in M and 
z to constrain Q m do not yet exist. Since the Press-Schechter mass function already assumes that 
X-ray clusters are virialized, one may convert this mass function to a temperature function with 
the virial theorem (see §2); however, X-ray cluster temperature catalogs that span sufficiently 
broad ranges in T and z to strongly constrain Q m also do not yet exist (Viana h Liddle 1998; 
Blanchard, Bartlett, & Sadat 1998; however, see Henry 1997; Eke et al. 1998). However, several 
X-ray cluster luminosity catalogs span sufficiently broad ranges in L and z to strongly constrain 
£l m , and the number of such catalogs is growing. However, to fit the Press-Schechter mass function 
to such catalogs, one must invoke a luminosity-temperature (L-T) relation in addition to the 
virial theorem. Theoretically, a wide variety of L-T relations have been proposed (e.g., Kaiser 
1986; Evrard & Henry 1991; Kaiser 1991); consequently, the L-T relation should be determined 
empirically. Until recently, the L-T relations of temperature catalogs have suffered from much 



-3- 



scatter (e.g., Edge & Stewart 1991; David et al. 1993; Mushotzky & Scharf 1997); however, 
recently Markevitch (1998), Allen & Fabian (1998), and Arnaud & Evrard (1998) have published 
temperature catalogs with temperatures and luminosities that have either been corrected for, or 
avoided the effects of cooling flows; the result is a significant reduction of this scatter. Using 
Bayesian inference, Reichart, Castander, & Nichol (1998) have constrained the slope and the 
evolution of the empirical L-T relation for the luminosity range 10 44 5 erg s _1 <; <J 10 46 - 5 erg 
s~ 4 and the redshift range z <J 0.5 from the Markevitch (1998) and Allen & Fabian (1998) catalogs. 

This latter work may be the key to determining cosmological parameters with X-ray cluster 
catalogs: given a well-constrained L-T relation, the Press-Schechter mass function may be fitted 
to the growing number of independent, high-redshift, high-luminosity X-ray cluster catalogs with 
well-understood selection functions. Using the cooling flow corrected L-T relation of Reichart, 
Castander, & Nichol (1998), we do this to the EMSS and the ROSAT Brightest Cluster Sample 
(BCS) luminosity function of Ebeling et al. (1997) in §3. In §2, we model the X-ray cluster 
luminosity function; in §4, we draw conclusions and discuss future applications of this luminosity 
function to the Southern SHARC and the Bright SHARC. 



2. The Model 

2.1. The X-ray Cluster Luminosity Function 

Assuming that X-ray clusters correspond to virialized, dark matter halos, we model the 
comoving number density of X-ray clusters with the Press-Schechter mass function, which is given 
by (e.g., Lacey k Cole 1993) 
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where n c (M,z) is the comoving number density of X-ray clusters of mass M at redshift z, 5 c q{z) 
is the present linear theory overdensity of perturbations that collapsed and virialized at redshift 
z, (To(M) is the present linear theory variance of the mass density fluctuation power spectrum 
filtered on mass scale M, and p~o is the present mean mass density of the universe. In the case of 
zero cosmological constant, which we assume throughout this paper, the present overdensity is 
given by (Lacey & Cole 1993; Peebles 1980) 
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and 



We assume a scale-free mass density fluctuation power spectrum of power law index n, so the 
present variance is given by (e.g., Lacey & Cole 1993) 



3+n 

M\ — 6" 



.o(M)=^-j , (7) 

where as is the amplitude of the mass density fluctuation power spectrum over spheres of radius 
8/i _1 Mpc, and Mg is the mean mass within these spheres. 

We now convert equation (jl|) from a mass function to an appropriately defined luminosity 
function. Following the notation of Mathiesen & Evrard (1998), we begin by assuming that X-ray 
clusters' bolometric luminosities scale as power laws in mass and redshift: 

L w «Af(l + z) s . (8) 

As did Henry et al. (1992) in the case of the EMSS, we find that the fractions of this luminosity 
that fall into the EMSS band of 0.3 - 3.5 keV and the BCS band of 0.1 - 2.4 keV are well 
approximated by power laws in X-ray cluster temperature: 

Lx = fxT- p L boU (9) 

where fx = 0.989 ± 0.014 and (3 = 0.407 ± 0.008 for the representative temperature range of the 
EMSS (3 £ T £ 10 keV), and fx = 1-033 ± 0.012 and = 0.472 ± 0.008 for the representative 
temperature range of the BCS (1.5 <^T <^ 12 keV), where temperature is measured in keV. At 
lower temperatures, however, these approximations quickly fail. Equation @ is independent of 
redshift because X-ray cluster luminosities are measured in the source frame. The temperature 
dependence introduced by equation (||) is removed with the virial theorem: 

TocMI(1 + z). (10) 

Technically, this expression holds only when £l m = 1; however, we show in §4 that generalizing 



this expression has little effect upon our results. Together, equations (g), (||), and (|i0|) yield the 



following expression that relates an X-ray cluster's mass to its observed luminosity, Lx'- 

L x <xf x M*>- 2 -i(l + z) s -P. (11) 



- 5 - 



Substitution of equation ([Tl|) into equation (|l|) yields the following luminosity function: 
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and a and c depend upon as and the factor of proportionality of equation (11). Instead of trying 
to model this factor of proportionality and fitting to erg, we simply fit to such degenerate, or 
grouped, combinations of these parameters in this paper (see §2.2). 

Since luminosities are computed from measured fluxes and redshifts, Lx is a function of Hq 
and U m . With one exception, all dependences upon Hq can be grouped into the parameters a 
and c, and this exception is noted in §2.2. The EMSS and the BCS provide luminosities in their 
respective X-ray bands that have been computed for Hq = 50 km s -1 Mpc -1 and Q m = 1; we 
denote these luminosities by L\. The relationship between Lx and L\ is given by 



where 



L x =x(z)L 1 , (15) 



fW) 2 f F (d A (n m =i)) , Einstein) 



^l(^to) is luminosity distance, c2a(O to ) is angular diameter distance, and is the 

fraction of an X-ray cluster's flux that is detected in the 2'. 4 x 2'. 4 detect cell of the original 
EMSS.[] A complete description of this quantity can be found in Henry et al. (1992). Since ROSAT 
measures total fluxes, /^(^(^m)) = 1 here. So in the case of the BCS, the latter expression 
applies. However, the revised EMSS subsample that we fit to in §3.3 is a combination of 43 
Einstein luminosities and 18 ROSAT luminosities (see §3). Fortunately, 36 of the 43 Einstein 
clusters have redshifts of z < 0.33, and 6 of the remaining 7 clusters have redshifts of z < 0.47. 
At these redshifts, the ratio of fractional fluxes in the former expression for x(z) is within a few 
percent of unity for a wide range of values of £l m . Furthermore, luminosities are available 

for 7 of the 8 clusters that carry the majority of the weight in the fits of §3.3, and the remaining 
Einstein cluster is at a redshift of z = 0.259. Consequently, we also use the latter expression for 
x(z) in the case of the EMSS. Besides, we show in §3 that the sensitivity of x(z) to Q m plays only 



5 The quantity /p(dyt(fim)) is also a function of X-ray cluster core radius, ao, which Henry et al. (1992) found to 
be ao ~ 0.25 Mpc, assuming that Q m ~ 1. Repeating their analysis for Q m ~ 0, we find that this result again holds; 
consequently, we adopt this value of ao throughout this paper. 
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a tertiary role in the determination of this parameter. Substitution of equation (15) into equation 
( |i~2| ) allows the luminosity function to be fitted to L\ data without loss of generality: 
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2.2. The Selection Function 

Let A(Li, z) be the area of the sky that an X-ray survey samples at redshift z as a function 
of luminosity L\. In the case of the EMSS, this quantity is given by (Avni & Bahcall 1980; Henry 
et al. 1992; Nichol et al. 1997) 

A(L 1 ,z) = A(F hm = F(L 1 ,z)), (20) 

where A(Fu m ) is the area of the sky that the EMSS surveyed below sensitivity limit Fu m (see 
Henry et al. 1992), 

, , _ f F {d A {z)) fe§ Li 

F(Li ' z) - k( Z ) ^awr { ) 

and k(z) is the k-correction from the observer frame to the source frame for a T = 6 keV X-ray 
cluster; the exact dependence of k(z) upon X-ray cluster temperature can be ignored for the 
representative temperature range of the EMSS. For the EMSS, we have computed A(Li,z) for 41 
values of L\ between 10 43,5 and 10 45,5 erg s _1 (0.3 - 3.5 keV) and for Q m = 0, 0.5, 1, and 1.5. For 
intermediate values of L\ and f2 m , we use linear interpolation between 43.5 < logLi < 45.5 and 
< Q m < 1.5. The cases of £l m = and 1 are plotted in Figure 1. The dependence of A(Li,z) 
upon Hq cannot be grouped into the parameters a and c, unlike all of the other Hq dependences in 
this analysis (§2.1). Instead of making Hq a free parameter, we fix Hq = 50 km s _1 Mpc -1 in this 
paper. However, if others wish to be more general, they need only consider the Hq dependence of 
the selection function. The case of the BCS is treated separately in §3.2. 

The total number of X-ray clusters observed between luminosity and redshift limits 
L[ < L\ < L u and zi < z < z u , i.e., the cumulative luminosity function, is given by 

N{L h L u -z h z u ) = r A(L 1 ,z)^^dL 1 dV(z), (22) 
JLi Jzi dLi 

where 

, T/ , , _ 4c 3 dz (Q m z + (n m - 2)({Q m z + I) 1 ! -I)) 2 

dV{z) ~ Hl^ m (i + zf (i + n rnZ f 
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is the comoving volume element. Hence, our model (luminosity function + selection function) 
consists of nine parameters: Hq, fx, P, p, s, a, c, n, and fi m . We have fixed the value of Hq, and 
we have tightly constrained the values of fx and (5 for the EMSS and the BCS (§2.1). In §3.1, 
we adopt values of p and s from the Bayesian inference analysis of the X-ray cluster L-T relation 
of Reichart, Castander, & Nichol (1998), and in §3.3, the normalization parameter, a, drops out 
of the Bayesian inference analysis of this paper. This leaves three parameters: m , n, and the 
grouped parameter c, which depends upon cr 8 , the proportionality constant of equation (|Tl|) , and 
Hq. We determine credible intervals for the values of these three parameters in §3.3. 



3. Bayesian Inference 



If properly applied, equation (|22| ) can be an effective probe of Q m . In this cumulative 
luminosity function, 5 c q(z), A(Li,z), x(z), and dV(z) depend upon U m . We now consider how 
sensitive each of these quantities is to f2 m . The comoving volume element dV(z), is approximately 
given by (equation (f23|)) 

^cJwr <*--°» , (24) 

where the indices apply in the redshift range 0.14 < z < 0.6. The luminosity conversion expression, 
x(z), is approximately given by (equation (|l~6|)) 

*(*) « 1 - l -^z, (25) 
and in equation (0), it is always raised to a power that is between ~ —1.5 and ~ 0.5. 



Consequently, dV{z) and x{z) contribute only weak dependences upon VL m to equation (f22|). The 
present overdensity is a stronger function of Q m (equation (^)): 

I 1.69(1 + z) (O m = 1) 



Since this expression appears to the second power in the exponential cutoff of equation (|17D, it 
contributes a significant dependence upon £l m to the cumulative luminosity function. For example, 
if the luminosity function is observed to cut off prematurely at higher redshifts, i.e., if there is a 
deficit of high-redshift, luminous X-ray clusters, then higher values of £l m are favored. However, 
if little or no evolution is manifest in the observed X-ray cluster luminosity function, particularly 
above L^|, then lower values of Q m are favored. 

The surveyed area, A(Li,z), contributes a different type of dependence upon Sl m to the 
cumulative luminosity function. In the case of the EMSS (see Figure 1), this dependence 



6 In this paper, refers very generally to those luminosities at which the luminosity function, modeled by equation 
appears to roll over from a power law to an exponential cutoff. 
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is negligible at low luminosities and redshifts. However, at luminosities ;> L+, A(L±,z) is a 
non-negligible, increasing function of Q, m at sufficiently high redshifts. For example, in the case 
of a L\ = 10 45 erg s _1 , z = 0.8 EMSS cluster, A(Li,z) is roughly twice as large in an Q m = 1 
universe than it is in an Q m = universe. Although this effect is suppressed by the exponential 
cutoff of equation ( |T7| ) above L*, this effect is amplified about L* by the fact that the luminosity 
function itself is a non-negligible, increasing function of Q, m at luminosities <^ L* at these high 
redshifts (see §3.3). Consequently, we find that an overabundance of high-redshift, ~ L± EMSS 
clusters favors higher values of £l m and not lower values of this parameter as is generally thought. 
We return to this idea in §3.3. 

Consider first the case of X-ray cluster luminosity data that lies within a narrow redshift band 
of effective redshift z e fj. Then, up to a factor of A{L\, z e ff)dV(z e ff)/dz, the integrand of equation 
(|22| ) (equation (|i7|)) is simply a power law in luminosity with an O m -dependent exponential cutoff. 
This exponential cutoff is a function of the parameters fx, P, P, n, and g e ff = g(z e ff), which 
itself is a function of O m (see below). We have already constrained the values of fx and ft (§2.1), 
and we adopt the Reichart, Castander, & Nichol (1998) value of p, as well as that of s, in §3.1. 
However, there are too few high-luminosity X-ray clusters to simultaneously constrain n and g e ff- 



Fortunately, n is also constrained by the low-luminosity, power-law limit of equation (|17|) for which 
data is more plentiful. Consequently, by fitting this luminosity function to data of this type, n 
and g e ff can be jointly constrained. 



By equation (|19|) , the parameter g e ff is a function of z e ff and the parameters fx, ft, P, s, 
n, c, and £l m . The effective redshift is a given and the parameters fx, ft, P, s, and n can be 
constrained as described above. However, since z ~ z e ff, a constant, the parameters c and Sl m are 
degenerate; consequently, Sl m can only be constrained if the value of c is otherwise known, i.e., if 
the values of a%, the factor of proportionality of equation (|ll]), and Hq are otherwise known (§2.2). 
Even in the event that X-ray cluster mass data is used instead of luminosity data, any fitted value 
of Q m will still depend strongly upon the assumed value of as , as well as upon the assumed value 
of n, since mass data is not yet plentiful enough for the Press-Schechter mass function to constrain 
these parameters. 

However, now consider X-ray cluster luminosity data that spans a breadth of redshifts. 
Instead of constraining the single parameter g e ff, one instead constrains a distribution of such 
parameters with redshift, i.e., g(z). The normalization of this distribution is c and its shape 
yields J7 m since the parameters fx, P, P, s, and n are otherwise constrained. Consequently, 
by fitting equation ( |22"|) to the EMSS, which spans a breadth of luminosities and redshifts, the 
parameters n, c, and O m can be jointly constrained regardless of the values of <7g and the factor of 
proportionality of equation (|ll|) (but not regardless of the value of Hq since A{L\, z) is a function 
of this parameter (§2.2)). We do this for the EMSS in §3.3. In §3.2, we better constrain the 
parameters n and c (actually, g e ff) with the local {z e ff ~ 0.1) luminosity function of the BCS. 
First however, we discuss the cooling flow corrected L-T relation of Reichart, Castander, &; Nichol 
(1998) and its implied values of p and s, in §3.1. 
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3.1. The L-T Relation 

The combination of equations (S) and (|l^) yields the L-T relation: 

3p 3p 

Lbol oc T~2~(l + z) s ~~ . (27) 

Reichart, Castander, & Nichol (1998) have constrained the slope and the evolution of equation ( p7|) 
using the cooling flow corrected X-ray cluster temperature catalogs of Markevitch (1998) and Allen 
& Fabian (1998), and Bayesian inference. For the luminosity ranges of the EMSS and the BCS, 
and the redshift range z <; 0.5, they find that p = 1.86~^q \q and s = (3.77 — 0.63il m )^°;22- However, 
when using the L-T relation to fit the Press-Schechter mass function to X-ray cluster luminosity 
catalogs (1) that are not cooling flow corrected, and (2) for which X-ray cluster photon count rates 
have been converted to fluxes and luminosities by assuming a T = 6 keV thermal bremsstrahlung 
spectrum, they find that one should use p = 1.77+0 it an( ^ s = (3-14 — 0.65fi m )lQgg. 

In the case of the EMSS, this latter case applies. In this paper, we adopt the best-fit values: 
p = 1.77 and s = 3.14 — 0.65$7 m . The parameter p is well-constrained, and in §3.3, we show that 
the uncertainty in the value of s does not significantly affect our results. In the case of the BCS, 
luminosities are not cooling flow corrected, but they are not determined by assuming T = 6 keV 
for each X-ray cluster; instead, luminosities are determined by additionally requiring that each 
X-ray cluster satisfy a L-T relation (Ebeling et al. 1997). From Figure 1 of Reichart, Castander, & 
Nichol (1998), it is apparent that the assumption of T = 6 keV more strongly affects the value of 
p than does the use of cooling flow corrected luminosities. Consequently, in this paper, we adopt 
p = 1.86 in the case of the BCS. We note, however, that such minor variations in this parameter 
do not significantly affect our results. Since the BCS is a local (z e ff ~ 0.1) catalog, the value of s 
is unimportant in this case (see §3.2). 



3.2. The ROSAT BCS 

The ROSAT BCS is a flux-limited sample of 199 bright X-ray clusters. A complete description 
of this sample and its selection criteria can be found in Ebeling et al. (1997). The redshift range of 
the BCS is z < 0.3; however, most of the BCS clusters have redshifts of z < 0.2, and the effective 
redshift (see §3) of the sample is z e ff ~ 0.1. Consequently, the BCS samples the X-ray cluster 
population of the local universe. Although such a sample may not have enough redshift leverage 
to adequately probe £l m , its large size makes it an excellent sample to constrain the parameters n 
and c (actually, g e ff, §3). These constraints can then be combined with the EMSS results of §3.3 
to better constrain J7 m . However, since the BCS is not yet publicly available, we settle here for 
a simplified analysis of the binned BCS luminosity function of Figure 1 of Ebeling et al. (1997), 
which we rep lot in Figure 2. 

Since the BCS spans a relatively narrow band of redshifts, we let f(z) = f(z e ff) = f e ff and 
g(z) = g(z e ff) = g e jf in equation (|i~7|). This approximation is reasonable, unless the value of 
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Vt m is high, in which case one expects a lower comoving number density of X-ray clusters in the 
highest-luminosity bins. This is because the highest-luminosity bins more strongly sample the 
highest-redshift BCS clusters than do the lower-luminosity bins. Given this approximation — that 
z = z e ff, a constant — equation (17) can only reproduce such a high-luminosity roll-over of the 
luminosity function by favoring an artificially high value of n. To safeguard against this potential 
bias, we first fit equation (|i~7| ) to all 12 luminosity bins, then to all but the highest-luminosity bin, 
then to all but the 2 highest-luminosity bins, etc., until the fitted values of n and c do not change 
appreciably from fit to fit. Also, since Ebeling et al. (1997) have already corrected this luminosity 
function for sample completeness, we set A(L±, z e ff) = 1. Here, we ignore the dependence that 
this quantity has upon Q m , which is a reasonable approximation since z e ff < 0.3. 

We find that only the highest-luminosity bin noticeably changes our results. When we fit 
equation ( |i~7|) to all 12 luminosity bins, we find that n = — 0.47^o^ and g e ff = 0.901q^ 4 , where 
we have assumed a flat prior probability distribution between —3 < n < and < g e ff < 10, 
and the likelihood function is given by e~ x / 2 .[| When we fit equation @ to all but the 
highest-luminosity bin, we find that n = — 1.83lg and g e ff = 1.20^QgQ. Ignoring additional 
high-luminosity bins does not appreciably change this result; consequently, in this paper, we use 
all but the highest-luminosity bin to determine a constraint from the BCS. We will explore what 
this highest-luminosity bin implies for the value of Q m in a later paper. In Figure 3, we plot the 
1, 2, and 3 a credible regions in the n — g e ff plane for both of the above fits. We do not plot 
credible regions in the n — c plane because, by equation (1£), the parameter c is a degenerate 
function of the parameters g e ff and O m (see §3); however, credible regions in the n — c plane are 
easily recovered, given a value of Sl m (see §3.3). In Figure 2, we also plot the best-fit luminosity 
functions of these fits. 



3.3. The EMSS 

As described in §3, the breadth of the luminosity and redshift ranges of the EMSS makes this 
catalog an ideal sample with which to probe f2 m . In Figure 4, we plot the L\ — z distribution of 
the revised EMSS of Nichol et al. (1997), as well as the z < 0.14 portion of the original EMSS 
X-ray cluster subsample. The solid curves are contours of constant sampled differential volume, 
i.e., A(L\, z)dV(z)/dz = constant. From left to right, these contours are equally spaced from 
zero (zero contour not shown). If the X-ray cluster luminosity function has not evolved over the 
redshift range of the EMSS, then at each luminosity, most of the observed X-ray clusters would be 
where most of the sampled differential volume is. This appears to be the case below L\ ~ 7 x 10 44 
erg s~ x (0.3 - 3.5 keV), which demonstrates agreement with the results of Collins et al. (1997), 
Nichol et al. (1997), Burke et al. (1997), Rosati et al. (1998), Jones et al. (1998), and Vikhlinin et 
al. (1998a), i.e., that the X-ray cluster luminosity function evolves only minimally, if at all, below 



7 See, e.g., Gregory & Loredo (1992) for an excellent discussion of Bayesian inference. 
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(see §4). However, as was originally found by Gioia et al. (1990a) (see also Henry et al. 1992), 
there appears to be a deficit of high-redshift X-ray clusters above L±. For example, given that 
six X-ray clusters were detected between 0.14 < z < 0.4 at L x > 7 x 10 44 erg s~ 4 (0.3 - 3.5 keV), 
one would expect ~ 9 X-ray clusters between 0.4 < z < 0.6 and ~ 21 X-ray clusters between 
0.4 < z < 0.9 at these luminosities if the luminosity function were not evolving, yet only 2 and 4 
clusters were detected in these redshift ranges, respectively. Hence, we find a factor of 4 — 5 fewer 
X-ray clusters at redshifts above z = 0.4 than below this redshift at these luminosities. This is 
only possible at the ~ 1 — 2% level with the no evolution model. This suggests that high values of 
£l m may be favored (§3); however, we offer alternative interpretations of this deficit in §4. 

By Bayes' theorem, the posterior probability distribution for Q m , n, and c, P(Q m ,n,c), is 
given by normalizing the product of the prior probability distribution and the likelihood function 
(e.g., Gregory & Loredo 1992). Here, we assume a flat prior probability distribution between 
< £l m < 1.5, — 3 < n < 0, and < c < 3. The likelihood function, £(Q m ,n, c), is given by (e.g., 
Cash 1979) 

Ntot 

£(p, m ,n,c)= Y[ P(Li ti ,Zi\0, m ,n,c), (28) 

8=1 

where P(Ln,Zi\£l m ,n,c) is the probability that the ith X-ray cluster fits our model, given values 
of Q m , n, and c. For our model (equation (|22j)), this probability is given by (e.g., Cash 1979) 

P (A , Z | Sim ,n,c)= — 4fcf> jm^pA . m 

N(Li,L u ;zi,z u ) dz dL x 

Since both d>Ti c {Jj\^ J clL\ and N(yLi^L u \, zi,z u ) are proportional to the normalization parameter 
a, our results are clearly independent of this parameter. 

In equation (|28|), Ntot is the total number of X-ray clusters in the same region of the 
L\ — z plane as that over which N(Li, L u ; zi, z u ) is defined. This region should be as broad as is 
reasonably possible and it need not be rectangular, as the simple integration limits of equation 
( |22| ) suggest. In this paper, we set L\ < L\ < L M = 10 45 ' 5 erg s -1 , where L\ is set by the limiting 
flux of the EMSS: F(L U z) = 1.33 x 10" 13 erg cm" 2 s" 1 (Henry et al. 1992); hence, L x is a function 
of redshift.^ Also in this paper, we set 0.14 = z\ < z < z u = 0.6. We exclude higher redshifts 
because (1) the L-T relation of Reichart, Castander, and Nichol (1998) is derived from z ^ 0.5 
X-ray clusters, so its accuracy should not be trusted at redshifts much in excess of this value, and 
(2) optical identification become more difficult at these high redshifts (§1). This excludes three ~ 
L* X-ray clusters, which reduces the total number of X-ray clusters in our sample from 64 to 61. 
Below, however, we repeat our analysis with z u = 0.9 to show that the exclusion of these three 
X-ray clusters does not unfairly bias our results. In Figure 4, we mark this region of the L\ — z 
plane with a dotted line; solid points are interior to this region. 



8 In computing Li(z), we set Q, m — 0. This is because if we were to use a higher value of fi m when defining this 
curve, the selection function, A(L\,z), would be undefined for luminosities and redshifts near this curve when the 
value of Q, m is lower than the curve-defining value (see equation (pl[)). 
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The posterior probability distribution for any two parameters, e.g., P(Q m ,n), or one 
parameter, e.g., P(Q m ), is given by marginalizing the posterior probability distribution for all 
three parameters, P(Q m ,n, c), over the other parameters (e.g., Gregory & Loredo 1992). 1, 2, 
and 3 a credible regions are determined by integrating the posterior probability distribution 
over the most probable region of its parameter space until 68.3%, 95.4% and 99.73% of this 
distribution has been integrated (e.g., Gregory & Loredo 1992). In Figure 5, we plot the 1, 2, 
and 3 a credible regions of the two-dimensional posterior probability distributions P(£l m ,n) (top 
panel) and P(n,c) (bottom panel). The 1 a credible intervals of the one-dimensional posterior 
probability distributions, i.e., P(£l m ), P(n), and P(c), are VL m = 0.96+Q'gg, n = —2.28+0 25) and 
c = 0.66^Q23- Hence, the EMSS favors high values of fl m and low values of n. Although these 
results by themselves are constraining — Sl m < 0.2 is ruled out at the 2.3 a credible level — by 
combining the likelihood function of the EMSS with the posterior probability distribution of the 
BCS (Figure 3), stronger constraints can be placed upon these parameters. 

We now determine the combined posterior probability distribution of the EMSS and the BCS. 
Let the likelihood function be that of the EMSS, as given by equation (|2~8|). However, instead of 
assuming a flat prior probability distribution for all three parameters, as we did above, only assume 
a flat prior probability distribution between < VL m < 1.5, and use the posterior probability 
distribution of the BCS — PBCsin-ideff) (Figure 3) — equation (|l9|) , and the effective redshift 
of the BCS, z e jf ~ 0.1, to determine the full prior probability distribution: PBcs{^m-,n,c). 
In Figure 6, we plot credible regions of the two-dimensional combined EMSS/BCS posterior 
probability distributions P(Q m ,n) and P(n,c). In the three left panels of Figure 8, we plot the 
one-dimensional combined EMSS/BCS posterior probability distributions. The dotted lines in this 
figure mark the 1, 2, and 3 a credible intervals. We find that VL m = 0.96^32, n = — 1. 861^34, and 
c = 0.54^qj2- Values of O m < 0.2 are ruled out at the 3.0 a credible level. 

To establish that the exclusion of the three z > 0.6, ~ L± EMSS clusters does not unfairly 
bias our results, we repeat this analysis for z u = 0.9. In Figure 7, we plot credible regions 
of the two-dimensional combined EMSS/BCS posterior probability distributions P($l m ,n) and 
P(n, c). In the three right panels of Figure 8, we plot the one-dimensional combined EMSS/BCS 

I r\ Q Q _l ri Q r-7 

posterior probability distributions. In this case, we find that Q. m = 0.93I ' 2 6) n = — 1-50_ ' 36 , and 
c = 0.48^qj2- Values of S7 m < 0.2 are ruled out at the 3.5 a credible level. Lower values of Q, m are 
not favored for three reasons. First of all, it is clear from Figure 3 that we have added a great deal 
of volume for which there are no EMSS clusters above L±. Secondly, for the region of the L\ — z 
plane occupied by these three clusters, the EMSS surveyed roughly twice as much area if Q m = 1 
than it did if Vt m = (§3). This is primarily because higher values of Q m imply lower luminosity 
distances, which imply higher fluxes for a given luminosity and, consequently, greater surveyed 
areas by equation (|20|). The third reason, also mentioned in §3, is that at luminosities ^ L*, the 
luminosity function itself is a non-negligible, increasing function of f2 m . At these luminosities, the 
redshift dependence of the luminosity function is dominated by the function f{z). This function is 



-13- 



approximately given by (equations (|TsJ) , (g), and (|l6|)) 




1 case. Although our 



(30) 



z u = 0.9 results are more constraining than are our z u = 0.6 results, for the reasons stated above, 
we feel less confident about these results than we do about our z u = 0.6 results. 

We now show that the uncertainty in the value of s does not significantly affect our results. 
In the three left panels of Figure 9, we plot the one-dimensional combined EMSS/BCS posterior 
probability distributions, once again for z u = 0.6, except that we have now used the —la value 
of s from Reichart, Castander, & Nichol (1998). In the three right panels of Figure 9, we plot the 
same distributions, but for the +1 a value of s from Reichart, Castander, & Nichol (1998). The 
effect of varying the value of s by ±1 a is a variation in the values of the fitted parameters by less 
than the extent of their ±1 a uncertainties. If we add these uncertainties in quadrature, we find 
that f2 m rj 1.0 ± 0.4, n ~ —1.9 ± 0.4, and c ~ 0.55^g '^5. One-dimensional credible intervals for all 
of the fitted values in this section are compiled in Table 1. 



In this paper, we have constructed from the Press-Schechter mass function and the empirical 
X-ray cluster L-T relation of Reichart, Castander, & Nichol (1998) an X-ray cluster luminosity 
function that can be applied to the growing number of independent, high-redshift, X-ray cluster 
luminosity catalogs to constrain cosmological parameters. In particular, we have incorporated 
the evolution of the L-T relation and all significant dependences upon Q, m of the luminosity and 
selection functions into our Bayesian inference analysis. For a fixed value of Hq, we have applied 
this luminosity function to broad subsets of the revised EMSS X-ray cluster subsample of Nichol 
et al. (1997) and to the ROSAT BCS luminosity function of Ebeling et al. (1997) to constrain 
tt m . For the 61 revised EMSS clusters between 0.14 < z < 0.6, we find that VL m = 0.96±g || 
and n = -1.86±£g; for all 64 revised EMSS clusters, between 0.14 < z < 0.9, we find that 

1 o no _i_n t x r T 

Q m = 0.93_ 26 an d n = ~ 1-50-036- These high values of Q m are the result of an apparent 
deficit of high-redshift, luminous X-ray clusters, which suggests that the X-ray cluster luminosity 
function has evolved above L+. 

Nichol et al. (1997) suggested that the statistical evidence for the evolution of the EMSS 
luminosity function was only minimal. At first glance, this appears to be in contradiction to one 
of the conclusions of this paper. However, since Nichol et al. (1997) used a power-law luminosity 
function - which they did for purposes of comparison with the original EMSS result of Henry 
et al. (1992) - instead of a luminosity function that permits different degrees of evolution below 
and above L+, as we have done in this paper, their results are most directly applicable below L±: 
this is the luminosity range of the vast majority of the EMSS clusters, so it is by this luminosity 



4. 



Discussion & Conclusions 



-14- 



range that their results have most strongly been weighted. The fact that the X-ray cluster 
luminosity function does not evolve below L± has since been shown by Collins et al. (1997), Burke 
et al. (1997), Rosati et al. (1998), Jones et al. (1998), and Vikhlinin et al. (1998a). That the 
luminosity function appears to evolve above L+ is in agreement with the original EMSS findings of 
Gioia et al. (1990a), as well as the findings of Vikhlinin et al. (1998a,b) with the 160 deg 2 survey. 

The value of Q m that we find that this high-luminosity evolution in the EMSS corresponds to 
is consistent with the values found by Sadat, Blanchard, & Oukbir (1998) (O m = 0.85 ± 0.2) and 
Blanchard & Bartlett (1998) (Q m ps 1), based upon the work of Oukbir & Blanchard (1992,1997). 
Our value of Q m is somewhat consistent with the values found by Henry (1997) (f2 m = 0.50 ±0.14) 
and Eke et al. (1998) (Q m = 0.45 ± 0.2); however, Viana & Liddle (1998) have performed a 
more extensive error analysis upon a conservative subset of the data of these authors and find 
that O m ~ 0.75 with £l m > 0.3 at the 90% confidence level and £l m ~ 1 still viable. Blanchard, 
Bartlett, & Sadat (1998) find almost identical results (O m ~ 0.74, with 0.3 < VL m < 1.2 at the 95% 
confidence level) from these data. Finally, our value of O m is inconsistent with the values found 
by Bahcall, Fan, & Cen (1997) (O m = 0.3 ± 0.1), Fan, Bahcall, & Cen (1997) (VL m m 0.3 ± 0.1), 
and Bahcall & Fan (1998) (VL m = 0.2±°;f). 

Our value of n is consistent with the values found by Henry & Arnaud (1991) (n = — 1.7^'^) 
and Henry et al. (1992) (n = — 2.10^015), where these authors set f2 m = 1. Our value of n is 
also consistent with the value found by Eke et al. (1998) (n = — 1.69^^), where these authors 
included Q m as a free parameter. Our value of n is somewhat consistent with the value that 
Bahcall, Fan, & Cen (1997), Fan, Bahcall, & Cen (1997), and Bahcall & Fan (1998) adopted 
(n = -1.4). 

Taken as an ensemble, these results are perhaps discouraging in that they span the entire 
range of acceptable solutions: 0.2 ^ Q m <^ 1. This suggests that as yet unknown systematic effects 
may be plaguing some, if not all, of these results. We briefly identify seven areas where systematic 
effects could enter ours and similar analyses. (1) The first is the Press-Schechter mass function 
itself; however, numerical simulations (e.g., Eke et al. 1996; Bryan & Norman 1997; Borgani et 
al. 1998) consistently show that the Press-Schechter mass function is an adequate approximation. 
(2) The spherical collapse model of cluster formation (equations @ and (p6|)) may be inadequate. 
For example, numerical simulations by Governato et al. (1998) suggest that in equation (|26|), 
the expression 1.69(1 + z) may really be as low as ~ 1.6(1 + z) ' 9 . This suggests that use of the 
spherical collapse model may lead to underestimated values of ft m ; however, this is only a <; 10% 
effect. (3) Technically, equation ( |l0| ) only holds when Q, m = 1. Recently, Voit & Donahue (1998) 
derived a virial theorem that holds for all values of f2 m , and that allows for the fact that clusters 
grow gradually; their mass-temperature (M-T) relation reduces to equation ( |i~0| ) when Q m = 1. 
We find that M-T relations with functional forms that are similar to that of the M-T relation of 
Voit &; Donahue (1998) reduce our fitted value of fl m by <^ 10%; however, further investigation 
and use of this M-T relation is clearly needed. (4) Also on the subject of the M-T relation, 
care must be taken when fitting to X-ray cluster temperature catalogs: cooling flows lower the 



-15- 



measured temperature of most X-ray clusters, which should systematically affect values of Vt m that 
are determined in this way. (5) Based upon the cooling flow corrected X-ray cluster temperature 
catalogs of Markevitch (1998) and Allen & Fabian (1998), Reichart, Castander, & Nichol (1998) 
determined an empirical L-T relation between measured luminosities and cooling flow corrected 
temperatures that holds for z <; 0.5 and for luminosities that are typical of X-ray cluster catalogs 
(see §3.1); however, more cooling flow corrected X-ray cluster temperature measurements are 
needed to determine what, if any, exceptions exist to this L-T relation, and to extend it to higher 
redshifts. (6) The art of determining an X-ray cluster catalog's selection function is a constantly 
improving science; modern selection functions are determined via extensive numerical simulations. 
An alternative explanation to our high-f2 m result is that the EMSS, for whatever reasons, missed 
many high-redshift, high-luminosity X-ray clusters beyond what is accounted for by their selection 
function (see §3.3). However, given that the EMSS detected many high-redshift, low-luminosity 
X-ray clusters, this seems to be an unlikely scenario. (7) Finally, our cosmological model may be 
inadequate. We did not investigate the effects of a cosmological constant in this paper; however, 
many authors have demonstrated that the inclusion of a cosmological constant has little effect 
upon the determined value of J7 m (see §1). Also, the effects of quintessent and other exotic 
cosmologies have not yet been investigated in this context. 

Some of these potential sources of systematic error can be safeguarded against. For example, 
the spherical collapse model, the M-T relation, and the L-T relation all have proportionality 
factors that are potential sources of systematic error. However, as we have shown in §2.1, all of 
these factors, as well as the parameter ag, group together, giving us our parameter c. Since we 
fit for c, these factors cannot bias our result. However, inadequate functional forms for these 
relations, as well as for the other functions listed above, can bias ours and others' results. 

In addition to further theoretical and numerical development of this formalism, only the 
continued construction of X-ray cluster catalogs will act to further resolve these issues. Fortunately, 
the number of X-ray cluster luminosity catalogs is growing rapidly. One such catalog is the 
Southern Serendipitous High-Redshift Archival ROSAT Catalog (Southern SHARC) (Collins et 
al. 1997; Burke et al. 1997). The redshift and luminosity ranges of the Southern SHARC are 
z < 0.7 and L 1 < 3 x 10 44 erg s" 1 (0.5 - 2.0 keV). Although the Southern SHARC does not span 
the luminosity range of the EMSS, it will provide a good consistency check of our EMSS results. 
Our analysis of this catalog is underway. 

Two similar X-ray cluster catalogs that can serve a similar purpose are the ROSAT Deep 
Cluster Survey (RDCS) (Rosati et al. 1998) and the Wide Angle ROSAT Pointed Survey 
(WARPS) (Jones et al. 1998). The RDCS spans the redshift and luminosity ranges z < 0.8 and 
Li < 3 x 10 44 erg s" 1 (0.5 - 2.0 keV). The WARPS spans the redshift and luminosity ranges 
z < 0.7 and L 1 < 2 x 10 44 erg s _1 (0.5 - 2.0 keV). An analysis of the RDCS is also underway 
(Borgani et al. 1998). 

The 160 deg 2 survey (Vikhlinin et al. 1998b) and the Bright SHARC (Romer et al. 1998), 



-16- 



a high-luminosity extension of the Southern SHARC that is currently under construction, span 
redshift and luminosity ranges that rival those of the EMSS. Consequently, these catalogs will 
provide strong, independent checks of the EMSS results. 

Finally, local (z e jf ~ 0.1) X-ray cluster catalogs, such as the ROSAT BCS, are of great 
importance. Although these samples do not have the redshift leverage to constrain cosmological 
parameters, their large sizes make them excellent samples to better constrain the parameters n 
and c. Samples like the EMSS, the 160 deg 2 survey, and the Bright SHARC do not have sufficient 
luminosity leverage to strongly constrain these parameters, which leads to weaker constraints upon 
the cosmological parameters. However, a simultaneous analysis of a local X-ray cluster catalog — 
as opposed to a local X-ray cluster luminosity function as we have used in this paper — and any 
of these high-redshift, high-luminosity X-ray cluster catalogs could lead to significantly improved 
constraints upon all of these parameters. 
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and J. M. Quashnock for valuable discussions. We are also very grateful to our anonymous referee, 
whose input has greatly improved this paper. D. E. R. is especially grateful to Dr. and Mrs. 
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Table 1. Fitted Values of the X-ray Cluster Luminosity Function 



Catalog(s) z u s a Q m n c P(£l m > 0.2) 
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-i.83±g : ?g 








EMSS 


0.6 


3.14 


- 0.65fi m 


0.96«;« 


OO+0-36 
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0.66«;« 


2.3 
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EMSS+BCS 


0.6 


3.14 


- 0.65fi m 


0.961°;*! 


-i-86lg:S 


0.54lg;S 
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EMSS+BCS 


0.9 


3.14 


- 0.65fi m 


0.93±° f 6 


1 - OU _0.36 


0.48l°:» 
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EMSS+BCS 


0.6 


2.26 


- 0.61O m 


0.69±°;|? 


_i 74+0.42 


0.48l°;S 
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EMSS+BCS 


0.6 


3.79 


- 0.59fi m 


i -, 7 +0.33 
J^- 1 ' -0.28 


i qa+O-34 

i - ao -o.3o 


0-66±g:?| 


3.2 


a 



a Reichart, Castander, & Nichol 1998. 
9eff = l-20_ 60 . 
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Fig. 1. — The area A(L±,z) of the sky that the EMSS sampled at redshift z as a function of 
luminosity L\. The solid curve is for Q m = and the dotted curve is for tt m = 1. From left to 
right, the curves correspond to L x = 10 43 - 5 , 10 44 , 10 44 - 5 , 10 45 , and 10 45 - 5 erg s" 1 (0.3 - 3.5 keV) 
(see §2.2). 

Fig. 2.— The ROSAT BCS luminosity function of Ebeling et al. (1997). The solid line is the best 
fit of equation (18) to all 12 luminosity bins. The dotted line is the best fit of equation (18) to 
all but the highest-luminosity bin (see §3.2). The dashed line is the best-fit Schechter function of 
Ebeling et al. (1997). 

Fig. 3. — The 1, 2, and 3 a credible regions of the posterior probability distributions of the fit of 
equation (18) to all 12 luminosity bins (dotted lines) and all but the highest-luminosity bin (solid 
lines) of the ROSAT BCS luminosity function (see §3.2). 

Fig. 4. — The L\ — z distribution of the EMSS clusters. From left to right, the solid curves 
are increasing, equally spaced (from zero, zero contour not shown) contours of constant sampled 
differential volume (see §3.3). A deficit of high-redshift X-ray clusters is apparent above L\ ~ 
7 x 10 44 erg s _1 . The dotted curve is the more conservative of the two regions over which we fit 
equation (23) in §3.3. Points interior to this region are solid. The top panel is for Q m = and the 
bottom panel is for £l m = 1. 

Fig. 5. — The 1, 2, and 3 a credible regions of the marginalized posterior probability distributions 
P(fi m ,n) (top panel) and P(n,c) (bottom panel) of the fit of equation (23) to the 0.14 < z < 0.6 
revised EMSS clusters (see §3.3). 

Fig. 6. — The 1, 2, and 3 a credible regions of the marginalized posterior probability distributions 
P(Cl m ,n) (top panel) and P(n,c) (bottom panel) of the fit of equation (23) to the 0.14 < z < 0.6 
revised EMSS clusters and the ROSAT BCS luminosity function (see §3.3). 

Fig. 7. — The 1, 2, and 3 a credible regions of the marginalized posterior probability distributions 
P(fi m ,n) (top panel) and P(n,c) (bottom panel) of the fit of equation (23) to the 0.14 < z < 0.9 
revised EMSS clusters and the ROSAT BCS luminosity function (see §3.3). 

Fig. 8. — The marginalized posterior probability distributions P(Q m ), P(n), and P(c) of the fit 
of equation (23) to the 0.14 < z < z u revised EMSS clusters and the ROSAT BCS luminosity 
function. For the three left panels, z u = 0.6; for the three right panels, z u = 0.9. The dotted lines 
mark the 1, 2, and 3 a credible intervals (see §3.3). 

Fig. 9. — The marginalized posterior probability distributions P(fi m ), P(n), and P(c) of the fit 
of equation (23) to the 0.14 < z < 0.6 revised EMSS clusters and the ROSAT BCS luminosity 
function. For the three left panels, we use the —la value of s from Reichart, Castander, & Nichol 
(1998); for the three right panels, we use the +1 a value of s from Reichart, Castander, & Nichol 
(1998). The dotted lines mark the 1, 2, and 3 a credible intervals (see §3.3). 



-25- 




L, [0.3 - 3.5 keV] (10 44 erg s 1 ) 



- 26 - 




-27- 




- 28 - 




